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Precise measurements of the in-plane microwave surface impedance of high-quality single crystals 
of Bi2Sr2CaCu2 08+ x and Tl2Ba2Cu06+5 are used to probe the relaxation time of nodal quasipar- 
ticles in the d-wave superconducting state through a two-fluid analysis of the microwave conduc- 
tivity. While this analysis requires us to posit a form for the frequency-dependent quasiparticle 
conductivity, we clearly demonstrate that the extraction of the relaxation rate is quite insensitive 
to the assumed shape of the quasiparticle spectrum. The robustness of the analysis is rooted in the 
oscillator-strength sum rule and the fact that we simultaneously measure the real and imaginary 
parts of the conductivity. In both Bi2Sr2CaCu208+i and Tl2Ba2Cu06+<5 we infer a linear tem- 
perature dependence of the transport relaxation rate 1/r and a small but finite zero-temperature 
intercept. The linear temperature dependence of I/r is in accord with expectations for weak elas- 
tic scattering in an unconventional superconductor with line nodes and a small residual density of 
states. The same analysis reveals an onset of inelastic scattering at higher temperatures similar 
to that seen in the YBa2Cu306+a; superconductors. Finally we extrapolate the two-fluid model 
over a range of frequencies up to five times the measurement frequency, where the extrapolation 
predicts behaviour that is qualitatively similar to terahertz conductivity data on Bi2Sr2CaCu2 0g+a; 
thin films. While relaxation rates in Bi2Sr2CaCu2 0g+a; and Tl2Ba2CuOe+5 are substantially higher 
than in YBa2Cu30s+i there are qualitative similarities between all three materials, and the dif- 
ferences can likely be attributed to varying levels of static disorder. We therefore conclude that a 
universal picture of quasiparticle scattering in the cuprates is emerging. 
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I. INTRODUCTION 

A controversial issue in the field of high tempera- 
ture superconductivity is the nature of the intrinsic 
quasiparticle damping. Given the structural similari- 
ties of the different cuprate compounds and the simi- 
lar doping dependence of many of their physical prop- 
erties it is natural to expect the dominant mechanism 
of quasiparticle scattering to be universal. The exper- 
imental situation, however, is currently unclear. Early 
terahertai and microwave conductivity^ measurements 
on YBa2Cu30g+a: revealed a rapid drop in the quasipar- 
ticle relaxation rate upon entering the superconducting 
state. PhotoemissioniM and terahertai spectroscopies of 
Bi2Sr2CaCu208+a; provide a different picture, with a re- 
laxation rate for the nodal quasiparticles that is much 
larger than that of YP^CuaOg+a;, and is of the order 
of the thermal energy over the full temperature range of 
the superconducting state. Combined with the simple 
energy-temperature scaling of the single-particle lifetime 
observed in the photoemission experiments, this has led 
to proposals that proximity to a quantum critical point 
is responsible for the unusual quasiparticle dynamics in 
Bi2Sr2CaCu208+a;. 6 fn order to draw a more universal 
picture it is important to establish to what extent the 
intrinsic quasiparticle relaxation in the superconducting 



state of YBa2Cu306+3; differs from that in other impor- 
tant cuprate materials such as Bi2Sr2CaCu20g+ 2; and 
Tl2Ba2Cu06+,5. This might also clarify the differences 
in the transport relaxation rate measured by the mi- 
crowave conductivity and the single-particle scattering 
rate measured by photoemission experiments. Experi- 
mental probes sensitive to charge dynamics that can cross 
this divide and be applied to both YBa2Cu30g+a: and 
Bi2Sr2CaCu20g+ 2; are therefore of great importance to 
further progress in this area. 

Equally intriguing are the dynamical properties at 
the lowest temperatures, which should be dominated 
by nodal quasiparticles scattering from static disorder. 
The effect of impurities on the transport properties of 
a (i-wave superconductor is in general a complicated 
problem, requiring in most cases a self-consistent treat- 
ment of multiple scattering of quasiparticles from impu- 
rities. The standard theoretical model used for these 
calculations is the self-consistent i-matrix approxima- 
tion (SCTMA))^ 8 ^ which treats quasiparticle scattering 
within the framework of d-wave BCS superconductivity 
and Fermi liquid theory. In this model two simple regimes 
of scattering exist. The first is the weak-scattering or 
Born limit, in which the energy-dependent relaxation 
rate 1/t(E) is determined by the total phase space for 
recoil. As a result, \/t(E) should be proportional to the 
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density of states N(E) , which in a clean d-wave supercon- 
ductor is linear in energy. In that case the thermally aver- 
aged relaxation rate l/r(T) should have a linear temper- 
ature dependence. Recent broadband microwave spec- 
troscopy of Ortho-II-ordered YBa2Cu3 6.52 has for the 
first time revealed clear signatures of Born-like scatter- 
ing, including a linear temperature dependence of the av- 
erage relaxation rate and conductivity spectra with up- 
wards curvature down to the lowest frequencies. 10 The 
other regime in which the SCTMA makes simple pre- 
dictions is in the unitarity limit, in which the impu- 
rity has a quasiparticle bound state at the Fermi en- 
ergy that acts as a zero-energy scattering resonance, 
with 1/t(E) oc 1/E. Calculations have shown unitarity- 
limit scatterers to be strongly pair-breakingZi£ and have 
been strikingly confirmed by experiments in which the 
controlled addition of Zn impurities drives a crossover 
from linear to quadratic temperature dependence of the 
London penetration depth More recent measure- 

ments of the frequency-dependent microwave conductiv- 
ity of high purity Ortho-I YBa2Cu3 06.993 have refined 
the experimental situation and presented a puzzle — a 
temperature-independent relaxation rate below 20 K. 12 
This is compatible with neither the Born nor unitarity 
limits and has prompted extensions of the SCTMA the- 
ory to include the effect of 'order-parameter holes', lo- 
cal suppressions of the pairing amplitude at the impu- 
rity site that act as an additional scattering channel^ 
In the last year the SCTMA has been further extended 
to take into account small-angle scattering ii» This latest 
work is extremely promising, as the theory now starts 
from realistic models of the disorder potential, including 
that due to impurities located out of the Cu02 planes. 
Among its successes is the ability to model a range of ex- 
periments, including microwave and thermal conductiv- 
ity, angle-resolved photoemission and scanning tunneling 
spectroscopy, using just one set of parameters to describe 
the concentration, location and scattering phase shifts of 
impurities^ 

In this paper we use electrical transport experi- 
ments at microwave frequencies to probe the quasipar- 
ticle lifetime in near-optimally doped single crystals of 
Bi2Sr2CaCu208+2; and Tl2Ba2Cu06+a ■ Both the mi- 
crowave conductivity a = o~\ — \oi and the d.c. resistiv- 
ity are sensitive to the same fundamental physical prop- 
erty — the quasiparticle transport relaxation rate. How- 
ever, unlike d.c. resistivity experiments, measurements 
at high frequencies can probe electrical transport deep 
within the superconducting state. Since we simultane- 
ously measure both the real and imaginary parts of the 
microwave conductivity we are well positioned to ex- 
tract the temperature-dependent relaxation rate using 
a two-fluid analysis. This has been a somewhat spec- 
ulative endeavor in the past, with only limited exper- 
imental data available to support the assumption of a 
Lorentzian conductivity spectrum ^ifiiiLiS The two-fluid 
description of the conductivity we adopt here is based 
on a phenomenological conductivity spectrum cxi(w) that 



has been carefully validated through detailed spectro- 
scopic measurements on YBa2Cu306+ a; '^ i ^ Neverthe- 
less, we take care to show that the qualitative form 
of the temperature-dependent relaxation rate inferred 
from this analysis is quite insensitive to the assumed 
form of the conductivity spectrum. That this is so is 
primarily due to constraints imposed by the oscillator 
strength sum rule. In this paper we report new data 
on Bi 2 Sr2CaCu208+2; at temperatures down to 2 K, 
which we analyze along with previously published data 
on Bi 2 Sr 2 CaCu20 8+2; and Tl 2 Ba 2 Cu0 6+5 a2i2i The high 
resolution of the new experiments allows the scattering 
dynamics of the nodal quasiparticles to be revealed in 
great detail. In Bi 2 Sr2CaCu208+a; and Tl 2 Ba 2 Cu06+5 
the dominant sources of disorder are known to be located 
away from CuC>2 planes and should therefore act as weak, 
small- angle scatterers. We find that the transport re- 
laxation rates in Bi2Sr2CaCu20s+a; and Tl2Ba2Cu06+5 
have a similar temperature dependence, with 1/t(T) oc T 
at low temperature, as expected for a rf-wave supercon- 
ductor in the weak scattering limit. This is similar to ob- 
servations made in Ortho-II YBa2Cu30g.52r^ although 
the absolute relaxation rates are substantially higher in 
Bi 2 Sr2CaCu208+2; and Tl 2 Ba 2 Cu06. At higher tem- 
peratures 1/t(T) increases more strongly with tempera- 
ture, also reminiscent of transport in the YBa2Cu306+a; 
superconductors. 12 Our main conclusion is that while re- 
laxation rates in Bi2Si'2CaCu208+ :E and Tl2Ba2Cu06+5 
are higher than in YBa2Cu30g+a;, there are qualitative 
similarities, with the differences most likely due to vari- 
ations in the level of static disorder. 



II. EXPERIMENTAL TECHNIQUES 

Bi2Sr2CaCu208+a; crystals were grown in Leiden us- 
ing the travelling-solvent floating-zone (TSFZ) method 2 ^ 
and annealed in air at 500° C for 3 days to produce sam- 
ples with T c — 87 K. The samples have low levels of 
impurities and sharp superconducting transitions. They 
show a high degree of crystallinity, with x-ray rocking 
curves and magnetic-torque measurements indicating a 
mosaic spread less than 0.2°. In addition, the crys- 
tals were shown to be free of extended defects using 
magneto-optic Kerr microscopy. The crystal used in the 
microwave experiment was obtained by repeatedly cleav- 
ing a larger piece from a TSFZ boule to obtain a small 
0.5 mm x 0.5 mm x 0.02 mm platelet with shiny, mirror- 
like surfaces and sharp edges. In order to test and ex- 
tend the generality of the conclusions we draw from this 
new data we have carried out the same analysis on previ- 
ously published data from TSFZ Bi2Sr2CaCu208+:r crys- 
tals grown in Tsukuba by Kadowaki and Mochiku22i22i 
and flux-grown Tl2Ba2CuOg+5 crystals prepared in Cam- 
bridge by Mackenzie and Tyler ,21*24 Structural disorder 
in these materials has been characterised by electron- 
probe microanalysis (EPMA), which reveals cation non- 
stoichiometry at about the 10% level, an intrinsic con- 
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FIG. 1: Penetration depth data for the Leiden 

Bi2Sr2CaCu2 08+:E crystal show AA(T) oc T at low temper- 
ature. Inset: 38.6 GHz surface resistance data for the same 
crystal shows a sharp superconducting transition at T c = 87 K 
and a low residual surface resistance of 2 mfi. 

sequence of thermal-equilibrium fluctuations in compo- 
sition at the growth temperature. The EPMA measure- 
ments give average compositions of Tli.gsP^Cui.isOe+s 
for the Tl-cuprate, indicating Cu substitution onto the Tl 
sites, and Bi2.i6Sri.giCai.o3Cu208+a; for the Bi-cuprate. 
Since this cation disorder is located away from the Cu02 
planes the in-plane charge dynamics will only be weakly 
influenced by the impurity potentials, which should give 
rise to predominantly small-angle scattering^ 

Measurements of the surface impedance Z s = R s + iX s 
were made in circular cylindrical TEon Nb cavity res- 
onators, held at the helium bath temperature (1.2 K or 
4.2 K) and supporting TEon modes with Q factors in 
the range 2 x 10 6 to 3 x 10 7 . During these measurements 
the sample is positioned inside the cavity on a movable 
and independently-heated sapphire rod at an antinode of 
the microwave magnetic field, which is oriented perpen- 
dicular to the a-b plane of the sample. Although this ge- 
ometry involves substantial demagnetising fields for our 
platelet samples, it results in primarily in-plane screen- 
ing currents being induced, an important requirement for 
highly anisotropic materials like P^S^CaC^Os+z and 
Tl 2 Ba 2 Cu06+,5. The surface impedance of the sample 
is obtained from the measured cavity response using the 
cavity perturbation formula— 

A/ B (T) - 2iA/ (T) = T(R S + iAX s ), (1) 

where A/b(T) is the change in bandwidth of the TEon 
mode upon inserting the sample into the cavity, A/o(T) is 
the shift in resonant frequency upon warming the sample 



FIG. 2: (Color online) The 38.6 GHz microwave conductiv- 
ity of the Leiden Bi2Sr2CaCu208+a; crystal inferred from the 
surface impedance data in Fig. <ji (T) is scaled up by a fac- 
tor of 10 for clarity. A characteristic feature of all microwave 
conductivity data on Bi2Sr2CaCu2 08+a; and Tl2Ba2Cu06+i 
is the finite zero-temperature intercept of <ji(T). 



from Tbase to T, and T is an empirically-determined scale 
factor that depends only on the geometry of the sample 
and cavity. A key feature of our method is that R S (T) 
and AA S (T) w lu(1qA\(T) are measured at the same 
time and on the same sample, allowing both <J\{T) and 
"2 (T) to be determined at the measurement frequency in 
a model-independent way. 



III. EXPERIMENTAL RESULTS 

Figure □ shows R S (T) and AA(T) = A(T) - A(2 K) 
for the Leiden Bi 2 Sr2CaCu208+a; crystal at 38.6 GHz. 
At T c — 87 K Rs (T) has a sharp superconducting transi- 
tion, indicating good sample homogeneity, then decreases 
monotonically below T c to a low residual value of 2 mQ. 
No discernible change in R s was observed as the input 
power of the resonator was varied over more than two 
orders of magnitude implying an absence of extended de- 
fects, which are expected to give nonlinearities associated 
with weak superconducting links.. At low temperatures 
AA(T) has a strong, linear temperature dependence, a 
hallmark of d-wave pairing in the cuprates^ 7 - only ob- 
served in samples with low levels of strong-scattering 
defects 4i 

The complex microwave conductivity is obtained from 
R s and AX S using the local-limit expression o~\ — 102 = 
iw/Lto/^s assuming, as in earlier workfSSthat A(T = 0) = 
2100 A. A different choice of A(T = 0) does not modify 
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the qualitative features of <j(T), or of l/r(T) obtained 
from the two-fluid analysis. The conductivity data are 
shown in Fig. with <Xi(T) expanded by a factor of 
ten. While similar in form to previously published re- 
sults on Tsukuba-grown F^S^CaC^Os+z crystals, the 
new data are of higher resolution and extend to lower 
temperature. At low temperature <J2{T) oc 1/A 2 (T) has 
a strong linear temperature dependence, a consequence 
of line nodes in the energy gap. The data near T c are 
reminiscent of YE^CusOg+x, for which the superfluid 
density approaches zero with vertical slope due to 3D-XY 
critical fluctuations^ Fluctuations are also responsible 
for the sharp peak in <Ji{T) just below T c , 29 which is fol- 
lowed at lower temperatures by a broad peak near 20 K 
before decreasing slightly to a zero-temperature intercept 
of 4.5 x 10 6 ^- 1 m" 1 . 



IV. TWO-FLUID MODEL 

In this section we develop a two-fluid model in order to 
extract the thermally averaged quasiparticle relaxation 
rate from the microwave conductivity data. The com- 
plex conductivity of a superconductor has two contribu- 
tions that conduct in parallel: er s from the superfluid 
condensate; and a qp from quasiparticle excitations. The 
superfluid conductivity is dissipationless and takes the 
form 



Ti 6 / I 

s (u) = - 2 — - tt8(u) + — 
to* \ iuj 



(2) 



The zero-frequency delta function is required by causal- 
ity and represents the energy absorbed in accelerating 
the superfluid. The overall scale of a s is parameterized 
by n s /m* , the ratio of the effective density of supercon- 
ducting electrons n s to their effective mass to*. This 
quantity is temperature dependent and is closely related 
to the London penetration depth Al: 
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n B (T) = 

to* fio^XKT)' 



(3) 



The quasiparticle conductivity has a more complicated 
form than the superfluid term. It contains information on 
both the quasiparticle excitation spectrum and the quasi- 
particle charge dynamics. In a d-wave superconductor 
the conventional approach to calculating the quasipar- 
ticle conductivity is the self-consistent i-matrix approx- 
imation, which accurately takes into account the multi- 
ple scattering of quasiparticles from impurities, a process 
that leads to pair-breaking and scattering resonancesA£*2i 
It has been shown that the SCTMA calculations of the 
conductivity, for any scattering strength, can be well ap- 
proximated by a simple, energy-averaged Drude form: 8 



2 f oo 

a(tu 7 T) = / N(E) 



df(E) \ 1 
dE ) \w + 1/t{E) 



dE, 
(4) 



where 1/t(E) is the energy-dependent electrical- 
transport relaxation rate, N(E) is the density of states 
of the (i-wave superconductor and f(E) is the Fermi 
function at temperature T. In a c?-wave superconduc- 
tor the relaxation rate takes a simple form in two im- 
portant limits: for weak scattering (the Born limit) 
1/t{E) oc E; while for strong scattering (the unitarity 
limit) 1/t(E) oc l/E, with a scattering resonance at 
zero energy. The energy dependence of the relaxation 
rate leaves a clear signature in the quasiparticle conduc- 
tivity spectrum, making microwave frequency measure- 
ments well suited to probing the quasiparticle dynamics. 

Several experiments have attempted to map out the de- 
tailed form of the quasiparticle conductivity spectrum in 
the cuprateSfi2*i2ii^ with the clearest to date being that 
of Turner et al. who used a novel broadband spectrom- 
eter to measure the microwave conductivity of Ortho-I 
YBa 2 Cu3 6.993 and Ortho-II YBa 2 Cu3 06.52 as continu- 
ous functions of frequency from 1 GHz to 21 GHz. These 
measurements revealed conductivity spectra that were 
broadly in accord with those expected from the SCTMA 
calculations, with upwards curvature down to the low- 
est frequencies It was shown by Turner et al. that 
the conductivity spectra were well described by a simple 
phenomenological form: 



1 + (ujt)v ' 



(5) 



For y = 2 this spectrum has the familiar Drude form of 
a single Lorentzian. As y is reduced below 2 the up- 
wards curvature extends to lower frequencies and the 
spectrum starts to correspond well to the expected form 
for a d-wave superconductor with energy-dependent scat- 
tering. In fits to experimental data in Rcf. [Hjit was found 
that for Ortho-II-ordered YBa2Cu3 O6.52 V ~ 1-45 and for 
fully oxygenated YBa2Cu3 06.993 y ~ 1.67. 

There are several caveats to note when using Eq.|SJ For 
y =/= 2 the spectrum applies only at positive frequencies. 
Also, y must be greater than one or the spectrum cannot 
be integrated, in violation of the oscillator-strength sum 
rule. Finally, at high frequencies a(uj) ~ 1/uj v , whereas 
physically we expect o{uj) to eventually go as l/u> 2 since 
in all cases the energy-dependent relaxation rate entering 
Eq. 01 will have an upper bound. Eq.|31can therefore only 
be approximate at the highest frequencies. Neverthe- 
less, it has provided extremely good fits to the measured 
spectra for YBa2Cu306+a; in the frequency range 1 GHz 
to 21 GHziiS In those experiments Eq. [SI also appeared 
to extrapolate very well outside the measured frequency 
range, as the conductivity spectral weight inferred from 
integrating fits of Eq. [SJto the YBa2Cu3C>6+:E data very 
accurately tracked the loss of superfluid spectral weight. 
Nonanalytic conductivity spectra have been proposed in 
several other cases to describe the unusual optical con- 
ductivity of cuprate and ruthenate materials in the nor- 
mal state^S^i*^ In that work, the form of the conduc- 



5 



0.8 



0.6 



0.4 



0.2 









20 40 60 

Temperature (K) 



80 



6 
5 




,'a 4 



20 40 60 

Temperature (K) 



80 



FIG. 3: (Color online) Superfluid fraction f a (T) from a two- 
fluid analysis (y = 2) of conductivity data from the Leiden 
Bi2Sr2CaCu2 08+a; crystal. Open triangles show f s (T) for 
(jq = 7.5 x 10 7 fi^m" 1 = ct 2 (T -» 0). This choice of ct 
leads to a unphysical downturn in f s (T) at low temperature. 
Solid triangles shows / S (T) for cto = 8.3 x 10 7 f2 _1 m _1 , a 
choice that allows for a small amount of residual conductivity 
spectral weight as T — > 0. / s is now monotonic in tempera- 
ture. 



FIG. 4: (Color online) Relaxation rate l/r(T) from a two- 
fluid analysis (y = 2) of the microwave conductivity of 
the Leiden L^S^CaC^Os+z crystal. Open triangles cor- 
respond to cto = 7.5 x 10 7 fi — 1 m _1 , solid triangles to <tq = 
8.3 x 10 7 Q _1 m _1 . Allowing for residual normal fluid in the 
choice of cto results in a small offset to l/r(T). The data re- 
veal a clear linear temperature dependence of l/r(T) at low 
temperature followed by a strong increase in scattering at high 
temperatures. 



tivity was 



Ode 



(l/T + iuj) a 



(6) 



Spectra of this sort suffer the same limitations as Eq. 
but in each case their use was motivated by good fits 
to experimental data. With these warnings in mind we 
will continue on and use Eq. in the analysis of our 
conductivity data. We will be careful to show, though, 
that our conclusions do not depend on the assumed shape 
of a(u>). 

With a phcnomenological spectrum that we can ex- 
trapolate outside the measured frequency range, the 
imaginary part of the quasiparticle conductivity can be 
obtained by a Kramers-Kronig transform: 



,qp 



Ode 



i + (u>T)y 



io-dcKK(wr) , 



where 

KK(wt) 



2lo 

-V 
k Jo 



1 + (LU'T) g UJ' 2 - LU 



(7) 



rdw'. (8) 



Eqs. [21 and[7| can be combined into a phenomenologi- 
cal model of the full microwave conductivity of a d-wave 
superconductor at finite frequencies, which we will later 



use to extract the temperature dependent relaxation rate 
of the quasiparticles. The model is 



ne 
rn* 



k 

io; 



f n ry sin(7r/y) 



1 



1 + {lot) v 



i KK(wt) 



Here the prefactor of the quasiparticle term is obtained 
by insisting that the frequency integral of Eq. [S] 



Ode 



o 1 + (wr)* 



dy = 



0"dc 



(10) 



be independent of y and equal to ^f n ne 2 /m*, as re- 
quired by the oscillator-strength sum rule. In the clean 
limit, %/t <§; 2 A, the sum rule ensures that the super- 
fluid and quasiparticle conductivity spectral weights sum 
to a temperature independent value, imposing the con- 
straint that / s + / n = 1- This has been shown to apply 
in the case of the YBa2Cu30g+x superconductors^ and 
should certainly be the case in our I^S^CaC^Os+a; and 
Tl2Ba2Cu06+,5 crystals, as the superconducting energy 
gap A is large and the samples clean. 

In the limit y = 2, Eq. [5] reverts to the standard Drudc 
two-fluid model, 



io; 1/r 



In 



(11) 
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which in the past has been used to carry out analy- 
ses broadly similar to those presented below, albeit in 
the absence of experimental evidence to support the use 
of a Lorentzian spectrumiSiiSiiLiSiS For y — 2 there 
are closed-form expressions for the relaxation rate and 
normal-fluid fraction^ 



1/t 

fn 



c - cr 2 {T) 

(l+^ 2 T 2 )fl 



a 2 (T) 



CO 



(12) 
(13) 



Here oq = ne 2 /m*uj represents the total conductivity 
spectral weight and in general may be different from 
(72 (T — ► 0) due to the pair-breaking effects of disorder. 
Making allowance for residual normal fluid at T = is 
an important part of the analysis that follows. 

For j / 2, closed-form expressions for the relaxation 
rate and normal fluid fraction do not exist, forcing us 
to invert Eq. [5] numerically. Setting y < 2 leads to con- 
ductivity spectra that are very different from the Drudc 
form, but we show in the next section that the quali- 
tative temperature dependence of 1/r is insensitive to 
the choice of y. This is a consequence of the oscillator 
strength sum rule and helps to make our analysis very 
robust. 



ANALYSIS AND DISCUSSION 
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FIG. 5: (Color online) Relaxation rates 1/t(T) of Lei- 
den Bi2Sr2CaCu2 08+:c, Tsukuba Bi2Sr2CaCu2 08+;c and 
Tl2Ba2Cu06+i , obtained from a two-fluid analysis using a 
Drude (y = 2) quasiparticle conductivity spectrum. The low 
temperature linear regime and the strong temperature depen- 
dence of 1/t(T) at higher temperatures are generic features 
of the three samples. 



We begin by applying the two-fluid analysis to the new 
conductivity data from the Leiden B^S^CaC^Os+z 
crystal. As mentioned above, the total conductivity spec- 
tral weight, parameterized by o~o, plays an important 
role in the analysis. We demonstrate this as follows. 
First, we assume that all conductivity spectral weight 
condenses into the superfluid as T — » and set cr by 
extrapolating cr 2 {T) to T = 0. In this case we get 
(j o = 7.5 x 10 7 fi _1 m _:L . On carrying out the two-fluid 
analysis for y = 2 we obtain the data shown as open 
triangles in the plot of superfluid fraction / s in Fig. [3] 
f s (T) increases monotonically with decreasing tempera- 
ture down to 10 K, then peaks and shows a rapid down- 
turn. The origin of this unphysical behaviour is the as- 
sumption that all spectral weight condenses, which can 
be relaxed by choosing a larger value of o~q, corresponding 
to an increase in the total conductivity spectral weight. 
There is some uncertainty in how to choose ooi and in 
the analyses that follow we set cto to the smallest value 
that results in a monotonic temperature dependence of 
/ S (T). Other choices of oo do not affect our conclusions 
and make no qualitative changes to the form of l/r(T). 
For the Leiden B^S^CaC^Os+z crystal we have chosen 
cr = 8.3 x 10 7 n- x m -1 . The effect on f s (T), which is 
given by the closed symbols in Fig. |3 is quite dramatic 
at low temperatures: / s is now monotonic in tempera- 
ture. An extrapolation of f s (T) to T = gives a resid- 
ual, uncondensed spectral weight equal to 11% of the 



total (see Table 1). All microwave and terahertz conduc- 
tivity measurements to date find residual conductivity 
spectral weight: in YBa2Cu3 0g+x, where the spectrum 
can be measured in great detail, the residual conductivity 
spectral weight ranges from 0.5% to 3%ii2*iS The higher 
value in Bi 2 Sr 2 CaCu208+a; is likely due to increased lev- 
els of static disorder. Theorists have had some success in 
modeling the uncondensed spectral weight using impu- 
rities with scattering phase shifts intermediate between 
the Born and unitarity limitsi^*^ 

The effect on l/r(T) of varying o~o is not nearly as 
pronounced as it is on f s (T). Fig. 0] shows l/r(T) in- 
ferred from the Leiden Bi 2 Sr2CaCu208+a; data using the 
same choices of oo as above. Increasing ao gives a small, 
approximately temperature independent shift of l/r(T). 
These data also reveal one of the main findings of this 
paper, the strong linear temperature dependence of re- 
laxation rate at low temperature, which is the behaviour 
expected for weak-limit impurity scattering in a d-wave 
superconductor with line nodes. For the physically rel- 
evant case in which not all spectral weight condenses, 
1/t(T) approaches a finite low-temperature intercept. 
This is also what would be expected in the weak scatter- 
ing limit, in which the energy dependent relaxation rate 
is determined by Fermi's Golden Rule and should reflect 
the presence of a residual density of states N(E — > 0) 
induced by pair breaking. Both the linear temperature 



7 



* Leiden Bi Sr CaCu 

2 2 2 8+x 

Tsukuba Bi Sr CaCu 

2 2 2 8+x 

° Tl Ba CuO 

2 2 6+5 



CO 

(XI 

O 



y = 1.45 




FIG. 6: (Color online) Relaxation rates l/r(T) of Lei- 
den Bi2Sr2CaCu2 08+a;, Tsukuba Bi2Sr2CaCu20s+i and 
Tl2Ba2Cu06+5, obtained from a two-fluid analysis using a 
non-Drude (y — 1.45) quasiparticle conductivity spectrum. 
l/r(T) is reduced by a factor of » 1.5 but is qualitatively the 
same as for the Drude-limit analysis in Fig. 



dependence of 1/r and a finite residual relaxation rate 
are also observed in Ortho-II YBa 2 Cu3 06.52«^ 

We now apply the same analysis to our previously pub- 
lished conductivity data on Tl 2 Ba 2 Cu0 6 +5 (T c = 78 K) 
at 35.9 GHz and Tsukuba-grown Bi 2 Sr2CaCu20 8+x 
(T c = 91 K) at 34.7 GHz. To demonstrate the insen- 
sitivity of the analysis to the detailed shape of the con- 
ductivity spectrum, we carry it out in the y = 2 (Drude) 
limit and for y = 1.45, corresponding to the spectral 
shape observed in Ortho-II YBa2Cu3 06.52- In all cases 
the two-fluid analysis reveals the presence of significant 
residual normal fluid at T = 0. We do not include plots 
of f s {T), which are similar to Fig. EI but summarize the 
results for the residual normal fluid fraction in Tabled 



y 


Leiden 
Bi2 Sr2CaCu2 Os+x 


Tsukuba 
B^S^CaC^Og+i 


Tl 2 Ba2Cu0 6+i 


1.45 


0.16 


0.15 


0.21 


2 


0.11 


0.10 


0.15 



TABLE I: Residual normal fluid fraction, 1 - f s (T -v 0). 

The temperature dependent relaxation rates for all 
three samples are plotted in Figs. El and [fj] for y = 2 and 
y = 1.45 respectively. A close-up of the low temperature 
region of the y — 2 data, including a comparison with 
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FIG. 7: (Color online) A close-up view of the low tem- 
perature region of Fig. with the Bi2Sr2CaCu208+a; and 
Tl2Ba2CuOe+i data plotted alongside relaxation rate data 
from microwave spectroscopy of Ortho-I YBa2Cu30e.993 and 
Ortho-II YBa2Cu 3 06.52 i 10 ' 12 ' 19 Ortho-II YBa 2 Cu 3 06.52 also 
has a linear temperature dependence of l/r(T) a low temper- 
ature. The rapid onset of inelastic scattering at higher tem- 
peratures is a common property of all the materials shown. 



YBa2Cu306+a; data, is given in Fig. H The first notice- 
able feature is that the plots for different y are qualita- 
tively the same, differing only by a temperature indepen- 
dent scale factor » 1.5. This insensitivity of the spec- 
tral width parameter, 1/t(T), to the detailed shape of 
the conductivity spectrum has its origin in the oscillator 
strength sum rule. The superfluid density, which is ob- 
tained in our experiment through the simultaneous mea- 
surement of <72 (T) , fixes the oscillator strength of the con- 
ductivity spectrum at each temperature. A measurement 
of the height of the conductivity spectrum, a\{oj — > 0), 
then allows a very good estimate of the width of the 
spectrum, l/r(T). The second aspect of the l/r(T) 
data we notice is that there are broadly two regimes of 
temperature-dependent behaviour. At low temperatures 
1/t(T) is approximately linear in temperature. As dis- 
cussed earlier, this is the expected behaviour of a d-wave 
superconductor in the weak-scattering limit. Unlike a 
normal metal, in which l/r(T) approaches a constant 
low-temperature limit due to elastic impurity scattering, 
clastic scattering in a c£-wave superconductor is tempera- 
ture dependent: the degree of elastic scattering is better 
characterized by the slope of l/r(T). At higher tem- 
peratures l/r(T) grows much faster than linearly with 
temperature. We take this to be the onset of inelastic 
scattering. A regime of weak elastic scattering followed 
by a rapid onset of inelastic scattering was first observed 
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FIG. 8: (Color online) Extrapolation of the 38.6 GHz mi- 
crowave conductivity data from the Leiden Bi2Sr2CaCu208+a; 
crystal to high frequencies, using the two-fluid parame- 
ters from the Drude (y = 2) analysis. The extrapolated 
data are similar in form to terahertz conductivity data on 
Bi2Sr2CaCu2 08+:E thin filmsj- which show clear signs of 
residual conductivity spectral weight in the zero-temperature 
limit. 



in Ortho-I YBa2Cu3 06.993>^ In that context Walker and 
Smith pointed out the importance of Umklapp processes 
in relaxing electrical currents, showing that in a e?-wave 
superconductor such processes become gapped in the su- 
perconducting state, leading to an activated exponential 
temperature dependence of the relaxation ratei 3 ^ The 
distinction between Umklapp and normal processes has 
since been included in calculations of quasiparticles scat- 
tering from antiferromagnetic spin fluctuations to obtain 
good agreement with thermal and microwave conductiv- 
ity dataii£*2£ The third feature of the data is the similar- 
ity of the magnitude of l/r(T) in the three materials, be- 
ing nearly identical in the two Bi2Sr2CaCu208+ 2 ; samples 
and only 30% smaller in Tl 2 Ba2Cu06+,5- This suggests 
an intrinsic origin of the scattering, and it seems rea- 
sonable to look to cation nonstoichiometry as the domi- 
nant source of disorder. In both Bi 2 Sr2CaCu208+a; and 
Tl2Ba2Cu06+,5, deviations from cation stoichiometry at 
the 10% level are a result of thermal equilibrium fluctua- 
tions in composition at the crystal growth temperature. 
(Cation nonstoichiometry in YBa2Cu306+a;, by contrast, 
is at the 10~ 4 level, due to the high chemical stability of 
the YBa2Cu306+a; structure.) One consequence of the 
cation nonstoichiometry is the creation of disorder lo- 
cated off the Cu02 planes. This smooths the disorder 
potentials experienced by electron quasiparticles, which 
then give rise to predominantly small angle scattering. 



FIG. 9: (Color online) Extrapolation of the 38.6 GHz con- 
ductivity data to high frequencies using two-fluid parameters 
from the y = 1.45 analysis. The qualitative form of the high 
frequency extrapolation is insensitive to assumptions about 
the detailed shaped of the quasiparticle conductivity spec- 
trum, making this a robust analysis. 



Recent work to include the effect of small- angle scatter- 
ing in the SCTMA theory has been very successful in 
describing experimental data and represents a major ad- 
vance in our understanding of the transport properties of 
the cupratesi 14 ! 38 

We finish the analysis and discussion by making a con- 
nection with terahertz experiments on Bi2Sr2CaCu208+a; 
from the Berkeley group£ The terahertz measurements 
span a wide enough frequency range that within the 
superconducting state they capture most of the low- 
frequency conductivity spectrum. These experiments 
were carried out on Bi2Sr2CaCu208+:E thin films, so the 
samples are in principle quite different from the single 
crystals used in the microwave experiments. Neverthe- 
less, if a universal picture of quasiparticle scattering dy- 
namics in the cuprates is to emerge, it should be pos- 
sible to understand both the microwave and terahertz 
data within a common framework of intrinsic inelastic 
scattering and sample-dependent disorder. To motivate 
that comparison we use the two-fluid parameters for the 
Leiden Bi2Sr2CaCu208+a: crystal to extrapolate the mi- 
crowave conductivity o~i{u), T) into a range of frequencies 
comparable to those of the terahertz experiment. The ex- 
trapolated microwave conductivity is plotted in Figs. |S| 
and El for y = 2 and y = 1.45 respectively. The frequen- 
cies used for the extrapolation are / = 50, 75, 100, 150 
and 200 GHz. These are a factor of 4 smaller than the fre- 
quencies used in the terahertz experiment, / = 0.2, 0.3, 
0.4, 0.6 and 0.8 THz, a choice that reflects the lower level 
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of scattering in the single crystals compared to the thin 
films. The extrapolated microwave data are qualitatively 
similar to Fig. 1 of Ref. 5, with a large residual conduc- 
tivity (Ji (T — > 0) at low frequencies that becomes small 
at the higher frequencies. The ability of Nunner et al. 
to fit our Bi2Si'2CaCu208+a; microwave data within the 
SCTMA framework using a combination of dilute strong 
scatterers and a high density of weak-scattering extended 
defects! 4 " suggests that a similar calculation might ac- 
count for the terahertz data. 



VI. CONCLUSIONS 

We have carried out a two-fluid analysis of high resolu- 
tion microwave conductivity data on I^S^CaC^Os+x 
and Tl2Ba2Cu06+a single crystals. We have carefully 
tested the sensitivity of our analysis to the assumed form 
of the quasiparticle conductivity spectrum, and find the 
procedure to be robust. In all the F^Si^CaC^Os+i and 
Tl2Ba2Cu06+<5 samples under investigation we observe a 
rapid collapse of the transport relaxation rate upon en- 
tering the superconducting state, in apparent contradic- 
tion of photoemission and terahertz spectroscopies!^^ 
We also report the observation of a linear temperature 
dependence of the relaxation rate, theoretically predicted 
for d-wave superconductors in the weak-scattering limit. 
1/t(T) increases more rapidly at higher temperatures in 
a way that is reminiscent of the strong onset of inelas- 
tic scattering in YBa2Cu306+a;. On extrapolating the 



two-fluid model to frequencies up to five times the mea- 
surement frequency the extrapolation predicts qualita- 
tively similar behavior to terahertz conductivity data on 
Bi2Sr2CaCu208+i thin films£ Our main conclusions are 
as follows. Relaxation rates in Bi2Sr2CaCu208+:r and 
Tl2Ba2Cu06+a single crystals are substantially higher 
than in YBa2Cu306+2;, although not as high as in 
Bi2Sr2CaCu208+a: thin films. The form of the scatter- 
ing, in particular the linear temperature dependence of 
1/t, suggests that the dominant elastic scatterers act in 
the Born limit — these can likely be attributed to cation 
nonstoichiometry that creates defects located away from 
the Cu02 planes. With the linear, low temperature term 
in 1/t(T) attributed to elastic scatterers, the remain- 
ing temperature dependence of 1/t in Bi2Sr2CaCu208+a; 
and Tl2Ba2Cu06+a bears a strong similarity to that of 
YBa2Cu306+£c, pointing toward a universal mechanism 
for inelastic scattering in the cuprates. 
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